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■ i. Summary of Proposed Work 

This is a technical report summarizing our progress in our program of contributing high 
quality Earth-based occultation observations to NASA's Planetary Data System (PDS) 
Rings Node. During our first year of funding, we selected five data sets for eventual inclu- 
sion in the PDS Rings Node.These were Uranus occultation observations obtained by the 
PI and co-workers from the IRTF of event stars U34 (26 April 1986), U1052 (5 May 
1988), U65 (21 June 1990), U7872 (25 June 1991), and U7808 (28 June 1991). In our 
original proposal, we described four tasks: 

* Data sets to a common format 

* Documentation of the occultation observations and associated calibrations 

• Calculation of the occultation geometry for each event 

• Establish prototype PDS templates 

As discussed in our renewal proposal, submitted 8 June 1993, we have completed the 
first three tasks, and are working on the fourth. As an indication of our progress to date, in 
this report we provide information about each of the data sets, their formats, the documen- 
tation, and the method used for reconstructing the occultation geometry. 

II. Occultation Observations 

The five data sets described below (all unpublished) were taken under a wide range of 
observing conditions and the occulted stars differed considerably in brightness. We have 
intentionally included both some of the best recent stellar occultation observations, as well 
as some that are marginally useful. In this way, we hope to identify problems early and to 
establish a flexible data storage format for the PDS that can handle both high- and low- 
SNR data. The main difficulties with noisy data concern proper normalization of the signal 
in the midst of strongly varying background signal, and unambiguous identification of real 
ring features. Even high SNR data can be problematical if the detector introduces spurious 
high frequency noise into the data, or if the time constant of the IR system is so long that 
ring features are significantly smeared out. In our descriptions, we will selectively highlight 
potential problems in preparing each of the five data sets for eventual submission to the 
PDS Rings Node. 

U34 (26 April 1986) 

These observations were obtained from the IRTF by Marc Buie. They provide a good 
example of the many steps involved in putting data into a common format and documenting 
the calibrations. The data were originally in GRABBER format, the native format for the 
IRTF high speed photometer. They were converted by Buie into FITS format, and pro- 
vided to the PI on 9-track tape. We then converted the FITS files to our own internal stan- 
dard data formats (which we call ts - rime series and hsp - high speed photometry) for use 
with our in-house plotting and model fitting software. All of the data sets in this report have 
been converted to ts or hsp format, and once we have worked out a final prototype PDS 
data template with Phil Nicholson and Mark Showalter, we will write a routine that will 
convert both of these data formats to the standard PDS format. 

In addition to converting the data to a standard format, additional calibrations needed to 
be performed. For example, because the oscillator driving the data system at the IRTF was 
not perfect, timing signals recorded as part of the data stream had to be used to calibrate the 
actual data interval per point. Nominally, this was 2 msec/pt, but because of clock drift, we 
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determined the more accurate rate to be 1.9999702 msec/pt. As another example, micro- 
phonic pickup by the IR detector introduced considerable ripple and resonance noise into 
the raw data signal, and this needed to be filtered out by smoothing and by low-pass filter- 
ing. This data set has a marginal SNR for ring analysis: it is useful for ring profile mod- 
elling and ring orbit model fitting, but only after considerable pre-processing. 

A compressed display of the original data set is shown in Fig. 1. There is considerable 
baseline drift, because these (like all of the data sets described here) were observed in the 
DC mode (rather than in chopping mode), and thus long-term background variations are 
not suppressed. Located in the figure are the ring and atmosphere occultations, as well as 
the calibrations used to determine the level of full stellar intensity. It is quite clear that an 
uninitiated user would have a great deal of difficulty interpreting this light curve without 
some documentation to identify specific features in the curves. The observing conditions 
during ingress (Fig. 2) were poor, and even in this higher resolution plot, it is difficult to 
identify individual ring features. Ring egress (Fig. 3) was more favorable, and even with 
wandering baselines, individual rings are 'distinct. An example of an atmospheric occulta- 
tion is shown in Fig. 4, which is accompanied by an adjacent calibration check. The PDS 
Rings Node will need to be able to document such atmospheric occultations when they are 
part of the ring occultation observations as well. 

The data processing necessary to produce clear ring profiles includes multi-notch and 
low-pass filtering, as is evident from the power spectrum of a section of the data set shown 
in Fig. 5. Note the pickup of both 60 Hz (and harmonics) as well as other strong peaks. 
With only modest filtering, some ingress features can be identified in the noise of Fig. 2, 
such as the y ring profile fitted with a square well diffraction model in Fig. 6. The egress e 
ring profile (Fig. 7) has high SNR but clearly shows residual high-frequency noise that 
has not been sufficiently filtered out. An important decision to be made when preparing 
these data sets for final archive is to what extent the observer/contributor should pre-pro- 
cess the observations, as opposed to simply supplying the raw observations with text doc- 
umentation of known problems. 

As an example of the data file formats we use internally, a generic template ts data file 
header is shown in Fig. 8, and the actual data header for U34 is shown in Fig. 9. In the 
final submitted version, we propose also to include a table giving a chronological identifi- 
cation of significant features in the data set such as bias level changes, sky checks, identifi- 
able ring features, regions of cloudiness, and so forth. The final format for such a catalog 
will depend somewhat on the PDS occultation data format, but it will be an essential part of 
the documentation. 

Using all of the ring features identifiable in the data set, we have reconstructed the oc- 
cultation geometry and the location of the occultation chord across the rings and planet, as 
depicted in Fig. 10. Note that this is a nearly central occultation, but we found no evi- 
dence in the data of a Uranian central flash, which in any case would have been of very low 
amplitude for this chord because of the nearly spherical shape of the projected planet. From 
the solution for the occultation geometry, we can compute the equatorial plane radius 
probed by the occultation rays as a function of time, shown in Fig. 11 for this event. Note 
that, since many of the Uranian rings are inclined, this radius scale will not be adequate for 
precise measurements of ring features, and it will be important to decide how to distinguish 
between true ring plane radius and equatorial radius in the final PDS products. 
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U1052 (5 May 1988) 

This occultation was considerably higher SNR than U34, but as illustrated in Fig. 12, 
it was a grazing event, missing the planet and the inner four rings entirely (we estimate that 
it passed a mere 100 km away from the a ring at closest approach!). The full light curve is 
shown in Fig. 13, where individual ring features can be seen interspersed among calibra- 
tion checks, and where the background level changed substantially but smoothly. These 
observations were obtained by Marc Buie as well, and were translated from FITS into our 
internal ts format using the simple UNIX shell scripts and header files shown in Figs. 
14-16. The high quality of the data is illustrated in the 5 ring ingress and egress ring pro- 
files shown in Figs. 17-18, where the diffraction fringes are quite evident, as is the faint 
broad ring interior to the main part of the 5 ring. There is still some unremoved ripple in the 
data, as seen in these figures, but midtimes of the model fits are accurate to a few msec, 
making them extremely useful for ring orbit determination. The results of the full set of ring 
profile fits are given in the table labelled as Fig. 19. It would seem to be useful to include 
such results as part of the PDS package, although the results should also eventually be in 
print in a refereed journal. Finally, Fig. 20 shows the reconstructed equatorial plane ra- 
dius as a function of time. It will need to be decided if each datum should be tagged with an 
equatorial plane radius, or whether a simple polynomial function should be fitted, with the 
equatorial plane radius being computable with PDS-supplied algorithms. 

U65 (21 June 1990) 

The occultation of U65 occurred under marginal conditions and the star was relatively 
faint, making this one of the poorest quality data sets still worth analyzing for ring orbital 
studies. The observations were made by R. French and M. Clark from the IRTF. The re- 
constructed occultation geometry is shown in Fig. 21 - note that all the rings were oc- 
culted, but there was no planet occultation. The ingress light curve (Fig. 22) is quite 
noisy, and even in the egress curve (Fig. 23) it is challenging to identify a single convinc- 
ing ring feature. It is unlikely that this data set will ever provide useful constraints on the 
abundance of inter-ring material, or accurate optical depths of the known rings, but ring 
event times are still obtainable from the data, and hence the observations are still of some 
value. The headers for the ingress and egress data are shown in Figs. 24 and 25. To il- 
lustrate the difficulty of finding ring features in the ingress data, take a look at Fig. 26, 
which shows a ring profile fit to Ring 4 (at about 9 seconds from the start of the figure). 
The fitted ring profile occurs at very close to the predicted time for this ring, but no great 
weight should be placed on the resultant fit. This raises the issue of how to flag such 
marginal events - this one is probably real, but at the edge of plausibility. 

The situation is improved for egress, where the Ring 4 profile is easily detected (Fig. 
27). Note the residual presence of ripples in the data due to high frequency noise in the 
data. It should be mentioned that these observations were made using a treble boost circuit 
that amplifies the high frequencies that would otherwise be suppressed by the IR electron- 
ics. Thus, it is dangerous to assume that the noise is white noise, and formal errors from 
least squares fits to ring are almost certainly misleading. As further evidence that good ring 
profiles can be found in these data, Fig. 28 shows the egress profile of the a ring. Fi- 
nally, in Fig. 29, we show the computed equatorial plane radius as a function of time. 
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U7 872 (25 J une 1991) 

Serendipitously, two high SNR occultations were visible from the IRTF within three 
days of each other in June of 1991. R. French observed these events with N. Chanover 
and E. Mason. The reconstructed occultation geometry of U7872 is shown in Figs. 30 
and 31. The full occultation light curve is shown in compressed form in Fig. 32, where 
bias level changes, sky checks, and guiding problems produce a welter of dips and kinks in 
the curve. A complete set of ingress and egress ring profile fits was performed, with the re- 
sults shown in Figs. 33 and 34. Varying background levels are evident in a few of the 
ingress ring profiles (ex: (3, 5, and e), but the egress profiles are very stable and high SNR, 
and are among the best ever recorded. 

This is an example of a relatively clean data set that should require little additional work 
to prepare for final submission to the PDS. It is not hard to see that there is an inverse rela- 
tionship between the quality of the data and the amount of effort required to document the 
observations, and at some point it will need to be decided whether or not the effort is justi- 
fied for the marginal data sets. 

U7808(28 June 1991) 

The final data set we have processed as part of this work is the occultation of U7808, 
whose event geometry (Figs. 35 and 36) was very similar to the U7872 occultation just 
three days earlier. Because of the near-grazing incidence of the inner ring occultations, the 
ring plane radial velocity was very low, and the ring profiles were thus of extremely high 
SNR, particularly for the faint inner three rings (4, 5, and 6). The general character of the 
overall occultation light curve (Fig. 37) resembles that of U7872 - guiding problems dur- 
ing ingress were solved during egress, and the resultant ring profile fits (Figs. 38 and 
39) show some baseline variations on ingress (q and y) but very stable baselines for 
egress, and stunning diffraction fringes on the Ring 6 profile. This data set should also be 
relatively easy to put in final form for the PDS, once the file format has been established. 

ill. Future directions 

Clearly, the next step is to decide upon a template for ground-based one-dimensional 
photometric occultation observations for the PDS. Among the issues that make this chal- 
lenging are the desire to have all data in a common format along with the competing desire 
for simplicity and clarity. Occultation data from the HST, IR array occultation observa- 
tions, and multi-channel visual observations all have special requirements, and it will take 
some care to insure that the adopted formats are not cumbersome to use for data sets that 
don't satisfy a strict set of conditions. Because of these important questions, we have post- 
poned constructing templates until the fall of 1993, when R. French will be visiting Phil 
Nicholson at Cornell for several weeks prior to the annual DPS meeting. The plan will be 
to develop a template and to convert one of the data sets (probably U7808) into this format. 
Then, at the DPS meeting, we will present these proposals to the PDS Rings Node advi- 
sory group for comment and criticism. Based on the results of these sessions, we will re- 
vise the template as needed and then proceed during the coming year with formal prepara- 
tion of the data for submission to the PDS. 
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/usr/local/mit- 


SIMPLE = F / 
BITPIX = integer / 
CTYPE = string / 
NAXIS = integer / 
NAXISn “ integer / 
GROUPS = T / 
GCOUNT = integer / 
PCOUNT = integer / 

FILETYPE= ts / 
FLTPVERS= 1 / 
VERSDATE= ”1986/08/20" / 
VERSAUTH= jle, rgf , jak / 

CLABLn = string / 
CDSCRn = string / 
CUNITn = string / 
CDELTn = double / 
CRPIXn = double / 

CRVALn = double / 
CTIMEn - "hh :mm: ss * sss" / 
CZONEn = string / 
CDATEn = fT yyyy/mm/dd” / 


PLABLm 

= string 

/ 

PDSCRm 

= string 

/ 

PUNITm 

= string 

/ 

PZEROm 

= double 

/ 

PSCALm 

= double 

/ 


GLABLj 

= 

string 

/ 

GDSCRj 

= 

string 

/ 

GUNITj 


string 

/ 

GZEROj 

= 

double 

/ 

GSCALj 

— 

double 

/ 

INSTRj 

= 

string 

/ 

TLSCP j 

= 

string 

/ 

OBSVRj 

= 

string 

/ 

GSTTMj 

= 

"hh : mm: ss . sss" 

/ 

GZONEj 

= 

string 

/ 

GSTDTj 

= 

"yyyy/mm/ dd" 

/ 

GLGTHj 

« 

"hh :mm: ss . sss" 

/ 

GOFILj 

= 

string 

/ 

GLFILj 

— 

string 

/ 

LOGINk 

= 

string 

/ 

COMPk 

= 

string 

/ 

PTHNMk 

= 

string 

/ 

PRGRMk 

= 

string 

/ 

TIMEk 

- 

"hh:mm: ss . sss" 

/ 

DATEk 

= 

" y y y y / mm / dd " 

/ 


astron/data_doc/fits_templates/ts_l .ft 

T if standard FITS, F otherwise 
number of bits per data pixel 
declaration type for C (short, etc.) 
number of axes (1 for time series) 
number of data pixels on axis n (base 1) 

T means that groups are present 

# of groups (channels + other) 

# of parameters preceeding a group 

ts is time series file 
version number of the filetype 
date that this version was created 
author (s) of this file version 

label for the implicit axis n 

description of parameter m: more info than label 
units for the implicit axis n 
increment on implicit axis of 1 pixel 
index of reference pixel on axis n (zero base) 
-0.5 for ref. pixel makes CRVALn + CTIMEn the 
"start time" (i.e. when data were turned on) 
added value for axis n reference pixel in CUNITn 
base time for reference pixel on axis n 
time zone for CTIMEn 

base date for reference pixel on axis n 

for a time axis: Absolute Time of ref. pixel = 

CDATEn + CTIMEn + CRVALn (each in own units) 

label of axis for parameter m 

description of parameter m: more info than label 

units for parameter m 

zero for parameter m 

scale factor for parameter m 

value = scale_factor * pixel + zero 

label of axis for group j 

description of group j: more info than label 

units for group j 

zero for group j 

scale factor for group j 

value = scale_f actor * pixel + zero 

instrument used for group j 
telescope used for group j 
observer for group j 

start time of original observations for group j 
time zone for GSTTMj 

date for start time of original observations (j) 
length of original observation record 
original file containing "data" in group j 
last file containing "data" in group j 


user login id for creation of kth file 
computer used for creation of kth file 
pathname of working directory for kth file 
name of program that created the kth file 
time for creation of kth file 
date for creation of kth file 


/thalesf^ra^ 


SIMPLE = 
BITPIX = 
CTYPE 
NAXIS 
NAXIS1 = 
GROUPS - 
GCOUNT = 
P COUNT = 
FILETYPE= 
FLTPVERS= 
VERSDATE= 
VERSAUTH= 
CLABL1 = 
CDSCR1 = 
CUNITl = 
CDELT1 - 
CRPIX1 « 
CRVAL1 = 
CTIME1 = 
CZONE1 = 
CDATE1 = 
GLABL1 = 
GUNIT1 = 
GZEROl = 
GSCAL1 - 
GSTTM1 - 
GZONE1 • 
GSTDT1 * 
GLFIL1 = 
LOGIN1 - 
COMP1 
PTHNM1 = 
PRGRM1 = 
TIME1 
DATE1 
END 


F 

16 

short 

1 

4800000 / # pts in u34.all.l 
T 
1 
0 

ts 

1 

"1986/08/20" 
jle, rgf, jak 
time 

u34 irtf grabber 
seconds 
0.002 
-0.5 
0.0 

"14:00:27.688" 

UT 

"1987/02/26" 

2.2 microns 
counts per axis unit 
0.0 
1.0 

"1987/02/26" 

UT 

"14:00:27.688" / from Marc Buie's letter of 
TT /usr5/uranus/u34/ irtf /raw /u34 . all . 1 " 
rf rench 
thales 

"/usr5/uranus/u34/ irtf/ raw /grabber" 
raw2avg,pro (IDL) 

11:14 

"1991/11/22" 


13 March 1987 



BE 







/thales5/uranus/ul052/irtf7raw/unpack.sh 



# unpack the tape of 12 May 1988 data from Marc Buie 

# 

TAPE- / de v / n rmt 1 h 

MAXFILES=4 

i=0 

while [ $i -It $MAXFILES ] 
do 

i='expr $i + 1 ' 
outf ile=FILE$i .data 

echo Reading FILE$i and storing data only as $outfile 
# skip over the FITS header and swap bytes to form data files 
/bin/dd if =$TAPE bs=2880 skip=l conv=swab of=$outfile 
Is -1 $outfile 
done 

echo All Done! 
mt -f $TAPE rewoff 
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/thales5/uranus/u 1 052/irtf/raw/u 1052_raw2ts.sh 


# convert long data files from tape of May 1988 data from Marc Buie 

# to ts format. Note that the header file was edited by hand. 

# See unpack. sh for commands used to extract FITS data 

# 28 November 88 - Dick French ( rf rench@kepler) 

# 21 January 1989 - revised headers 

# 23 January 1989 - revised output file names 

# 25 January 1989 - change id number of star 

# 25 January 1989 - change id number of star 

# 11/21/91 - copied stuff to thales but ran as is 


inf ilel-FILEl .data 
outfilel=. . /ts/ul052irtf.25.ts 
headerl=FILEl . ts header 


infile2=FILE2 .data 
outfile2=. . /ts/ul052irtf . 05 . ts 
header2=FILE2 . ts_header 

inf ile3=FILE3. data 
outfile3=. ./ts/ul052irtf.01.ts 
header3-FILE3 . t s header 


inf ile4=FILE4 . data 

outf ile4=. . /ts/ul052irtf . 002. ts 

header4=FILE4 .ts_header 

echo -n Creating $outfilel . . . 
cat $headerl $infilel > $outfilel 
echo 

Is -1 $outfilel 

echo -n Creating $outfile2 . . . 
cat $header2 $infile2 > $outfile2 
echo 

Is -1 $outf ile2 

echo -n Creating $outfile3 . . . 
cat $header3 $infile3 > $outfile3 
echo 

Is -1 $outfile3 


echo -n Creating $outfile4 . . . 
cat $header4 $infile4 > $outfile4 
echo 

Is -1 $outfile4 


■pm 

SIMPLE = 
BITPIX = 
CTYPE 
NAXIS 
NAXIS1 = 
GROUPS = 
GCOUNT = 
P COUNT = 
FILETYPE= 
FLTPVERS= 
VERSDATE= 
VERSAUTH- 
CLABL1 = 
CDSCR1 *= 
CUNIT1 = 
CDELT1 - 
CRPIX1 = 
CRVAL1 - 
CTIME1 - 
CZONE1 - 
CDATE1 = 
GLABL1 = 
GUNITl = 
GZER01 = 
GSCAL1 = 
TLSCPl = 
OBSVR1 = 
GSTTM1 = 
GZONE1 = 
GSTDT1 = 
GLGTH1 « 
GOFIL1 = 
GLFIL1 = 
LOGIN1 = 
COMP1 
PTHNM1 - 
PRGRM1 = 
TIME1 
DATE1 
END 


mim 


/thales5/uranus/ul652/irtf/raw/FILE4.ts_header 



F / 

16 / 

short / 

1 / 

3600000 / 

T / 

1 / 

0 / 

ts / 

1 / 

”1986/08/20" / 

jle, rgf, jak / 

time / 

ul052 IRTF / 

seconds / 

0.0020 / 

-0.5 / 

0.0 / 

”13:00:00.016" / 

UT / 

”1988/05/12" / 

2.2 microns / 

counts per axis unit/ 

0.0 / 

1.0 / 

IRTF / 

Marc Buie / 

”13:00:00.016" / 

UT / 

"1988/05/12" / 

" 02 : 00 : 00 . 000 " / 

FILE4 , data / 

FILE4 . data / 

rfrench / 

kepler / 

"/usrl /uranus /u4 6 /mko/ raw" 
dd and unpack. sh / name 


"18:23:00" 

"1989/01/21" 


T if standard FITS, F otherwise 
number of bits per data pixel 
declaration type for C (short, etc.) 
number of axes (1 for time series) 
number of data pixels on axis n(base 1) 

T means that groups are present 

# of groups (channels + other) 

# of parameters preceeding a group 
ts is time series file 

version number of the filetype 
date that this version was created 
author (s) of this file version 
label for the implicit axis n 

description of parameter m: more info than label 

units for the implicit axis n 

increment on implicit axis of 1 pixel 

index of reference pixel on axis n (zero base) 

added value for axis n reference pixel in CUNITn 

base time for reference pixel on axis n 

time zone for CTIMEn 

base date for reference pixel on axis n 

label of axis for group j 

units for group j 

zero for group j 

scale factor for group j 

telescope used for group j 

observer for group j 

start time of original observations for group j 
time zone for GSTTMj 

date for start time of original observations (j) 
length of original observation record 
original file containing "data" in group j 
last file containing "data" in group j 
user login id for creation of kth file 
computer used for creation of kth file 
/ working directory for kth file 
of program that created the kth file 
time for creation of kth file 
date for creation of kth file 
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/thales5/uranus/u65/ts/u65_i.hdr 


SIMPLE = F 
BITPIX = 16 
CTYPE = short 
NAXIS = 1 
NAXIS1 = 1916000 
GROUPS = T 
GCOUNT = 1 
PCOUNT = 0 
FILETYPE= ts 
FLTPVERS= 1 

VERSDATE= "1986/08/20" 

VERSAUTH= jle, rgf, jak 

CLABL1 = time 

CDSCR1 = u65 ingress 

CUNIT1 = seconds 

CDELT1 = 0.002 

CRPIX1 = -0.5 

CRVAL1 =0.0 

CTIME1 = "10:00:00.299" 

CZONE1 = UT 

CDATE1 = "1990/06/20" 

GLABL1 = 2.2 microns 

GUNIT1 = counts per axis unit 

GZEROl =0.0 

GSCAL1 =1.0 

GSTTM1 = "1990/06/20" 

GZONE1 = UT 

GSTDT1 = "10:00:00.299" 

GLFIL1 = "/usrl/tmp/u65/u65i rt2f 3 . ts . hdr" 

LOGIN1 = mclark 

COMP1 = thales 

PTHNM1 = "/usrl/tmp/ u65/ts" 

PRGRM1 = grabber2dat . pro 

TIME1 = 15:35 

DATE1 = "1990/08/17" 

END 





SIMPLE - 
BITPIX = 
CTYPE 
NAXIS 
NAXIS1 = 
GROUPS - 
GCOUNT « 
P COUNT = 
FILE TYPE— 
FLTPVERS= 
VERSDATE= 
VERSAUTH= 
CLABL1 - 
CDSCR1 = 
CUNIT1 * 
CDELT1 » 
CRPIX1 = 
CRVAL1 = 
CTIME1 = 
CZONE1 « 
CDATE1 - 
GLABL1 = 
GUNIT1 * 
GZEROl = 
GSCAL1 - 
GSTTM1 * 
GZONE1 * 
GSTDT1 « 
GLFIL1 = 
LOGIN1 - 
COMP1 
PTHNM1 = 
PRGRM1 = 
TIME1 
DATE1 
END 
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